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We reinvestigate the momentum-resolved single-particle spectral function of the Tomonaga- 
Luttinger model. In particular, we focus on the role of the momentum-dependence of the two- 
particle interaction V{q). Usually, V{q) is assumed to be a constant and integrals are regularized 
in the ultraviolet ‘by hand’ employing an ad hoc procedure. As the momentum dependence of the 
interaction is irrelevant in the renormalization group sense this does not affect the universal low- 
energy properties of the model, e.g. exponents of power laws, if all energy scales are sent to zero. 
If, however, the momentum k is fixed away from the Fermi momentum kp, with \k — /cf| setting 
a nonvanishing energy scale, the details of V(q) start to matter. We provide strong evidence that 
any curvature of the two-particle interaction at small transferred momentum q destroys power-law 
scaling of the momentum resolved spectral function as a function of energy. Even for \k — kpl much 
smaller than the momentum space range of the interaction the spectral line shape depends on the 
details of V{q). The significance of our results for universality in the Luttinger liquid sense, for 
experiments on quasi one-dimensional metals, and for recent results on the spectral function of one¬ 
dimensional correlated systems taking effects of the curvature of the single-particle dispersion into 
account (‘nonlinear LL phenomenology’) is discussed. 


I. INTRODUCTION 

A. Luttinger liquid universality and the 
Tomonaga-Luttinger model 

It is well established that the Tomonaga-Luttinger 
model (TLM)PI^with linear single-particle dispersion and 
a two-particle interaction potential V{q) which is finite 
at zero momentum transfer q = 0 forms the infrared 
fixed point under renormalization group (RG) flow of 
a large class of gapless one-dimensional (Id) models of 
correlated fermions.l^This is the essence of the much cel¬ 
ebrated Luttinger liquid (LL) universality.l^M^ It implies 
that the low-temperature thermodynamic properties as 
well as the low-energy spectral functions of a model be¬ 
longing to the LL universality class are equivalent to the 
ones of the TLM. Understanding the low-energy physics 
of the latter is th us of crucial importance. Fortunately, 
using bosonizatiorP^ it is possible to derive exact and 
closed analytical expressions for thermodynamic observ¬ 
ables such as the specific heat or the compressibility as 
well as for space-time correlation functions of the TLM. 
From the latter spectral functions can be computed by 
Fourier transform. 

The bosonization expressions for correlation func¬ 
tions of the TLM depending on position x and time 
t generically contain integrals over momenta. Within 
constructive bosonization, which is based on operator 
identities,SI3 these are naturally regularized in the ul¬ 
traviolet by the momentum-space range gc > 0 of the 
two-particle potential V{q). For time-dependent corre¬ 
lation functions the momentum integrals cannot be per¬ 
formed even if a specific form of U(g) is assumed.How¬ 
ever, one can show that the momentum dependence of 
the interaction is RG irrelevant This is employed to 


justify the following procedure: In the final expressions 
for the space-time correlation functions V (g) is routinely 
replaced by a constant. As a consequence the momentum 
integrals become divergent in the ultraviolet. These di¬ 
vergences are re gularized (‘by hand’) applying an ad hoc 
procedureWe already now emphasize that this reg¬ 
ularization is not unique. After these steps the momen¬ 
tum integrals can be performed and integral-free expres¬ 
sions for space-time correlation functions are obtained. 
In field-theory inspired phenomenological bosonization 
procedure^ the momentum dependence of the interac¬ 
tion is often neglected from the outset (even in the 
Hamiltonian). In correlation functions this leads to the 
same ultraviolet divergences as described above requiring 
again a regularization ‘by hand’. Similarly, the purely 
fermionic approach to the single-particle Green function 
G^{x,t) of the TLM by Dzyaloshinskii and LarkirP^ re¬ 
quires an ad hoc ultraviolet regularization. 


We thus emphasize that the integral-free expressions 
for a variety of time-dependent correlation functions 
which can be found in the literature cannot be consid¬ 
ered as the corresponding exact correlation functions of 
the TLM. This is often acknowledged by stating that the 
integral-free expressions of the ad hoc procedure only 
agree to the exact ones at ‘asymptotically large space- 
time distances’; as discussed in Ref. |9] (see also below), 
even this is incorrect when considering the decay in the 
directions specified hy x = ±vt, with v being one of the 
eigenmode velocities at small momentum. 
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B. Spectral functions of the Tomonaga-Luttinger 
model — the fate of power laws 


We now focus on the two-point correlation function - 
the Green function - at temperature T = 0 from which 
the single-particle spectral function can be computed by 
Fourier transform. The spectral function is of partic¬ 
ular importance as it provides forthright access to cor¬ 
relation effects and can directly be related to photoe¬ 
mission spectra. We consider the momentum integrated 
function p^{oj) [p^(co’)], which is experimentally acces¬ 
sible in angular integrated [inverse] photoemission, as 
well as the momentum resolved spectral function p^{k^uj) 
cj)]. A measurement of the latter requires momen¬ 
tum resolution. It was shown that the universal low- 
energy power-law suppression of p^{uj) ^ \oj\^ for cj ^ T 
with a > 0 is unaffected by the above ad hoc procedure.^ 
Regardless of the details of V(q) the exponent a depends 
on the potential at vanishing momentum transfer R(0) 
(only), that is the constant interaction strength after the 
ad hoc procedure. Furthermore, the power-law nonana- 
lyticity of p^(kF,u;) ^ exactly at the Fermi mo¬ 

mentum k = kp remains the same; depending on the size 
of R(0) a divergence (a < 1) or suppression (a > 1) 
might occur. These findings are consistent with the RG 
irrelevance of the momentum dependence of V(q) as in 
both cases all energy scales, that is uj and, in the case of 
momentum resolved spectra, v^ik — kF) are sent to zeioW 
Here v-p denotes the Fermi velocity. 

The question we address here is whether or not any 
of the standard ad hoc procedures is legitimate when it 
comes to p^{k,uj) as a function of uj at fixed k — k^ 7 ^ 0. 
Employing these to compute p^{k,uj) of the spinless 
TLM, characteristic algebraic threshold nonanalyticities 
at :^Vc[k — k-p] with the charge velocity Vc of the (col¬ 
lective, bosonic) charge eigenmodes were found. In the 
model with spin additional algebraic nonanalyticities ap¬ 
pear at :^Vs[k — k^] with the velocity Vs of the spin 
modes.lIS®3l xhe corresponding exponents can be ex¬ 
pressed in terms of the (momentum independent) inter¬ 
action potential. For small interactions and k — kp < 0, 
p^{k,uj) shows power-law singularities at uo = Vc[k — kp] 
and UJ = Vs[k — /cp] instead of a single (Lorentzian) peak 
which would emerge in a Fermi liquid. This is one of the 
signatures of spin-charge separation regarded as a hall¬ 
mark of (spinful) LLs. In Ref.|T5]it was shown that these 
features can be found in the exact spectral function of 
the TLM if a box-like potential V{q) = V{0)Q{q‘^ — q^) 
is assumed, as long as \k — kp\ < q^ with the momentum 
transfer cutoff gc- We emphasize that using a box poten¬ 
tial is per se not equivalent to the ad hoc procedure, as it 
can e.g. be seen considering p^{k,uj) for \k — kp\ > qc^ 
For \k — kp\ < qc the box potential might, however, be 
viewed as a unique realization of the ad hoc procedure 
(see below). Clearly, a box potential is rather special 
and might not even be considered as particular physical. 
Thus further work for more generic V (g) is required. 

In Ref. [ 9 ] it was shown that the algebraic properties of 


the Green functions G^{x,t) in the space-time plane are 
significantly affected by the ad hoc procedures. It was 
proven that the exponent of the asymptotic decay of the 
Green function in the distinguished directions x = 
is not only set by V (0) but in addition by a measure of the 
flatness of the potential at g ^ 0 ; a result which cannot 
be obtained within any ad hoc procedure. Based on this 
and the crucial insight that the decay of G^{x,t) in the 
distinguished directions plays a central role in obtaining 
the power-law nonanalyticities in p^(k,uj), the question 
was posed if for generic V{q), p^{k,uj) is characterized 
by the ‘thought to be universal’ power laws of the ad hoc 
procedure. However, Ref. [9] lacks a definite answer. 


The above mentioned box potential is ‘infinitely flat’ 
at g -> 0 and thus ‘nongeneric’. The asymptotics of 
G^{x,t) for this and the ad hoc procedures agree and 
consequently also the features of p^{k, uj) at uj ^ '^c/s[k — 
kp]. 


It is crucial to realize that a dependence of the line 
shape of p^{k,uj) on the details of the interaction away 
from g = 0 including such fundamental issues as the 
presence or absence of algebraic nonanalyticities does 
not contradict the RG irrelevance of the momentum de¬ 
pendence of V{q) in the TLM. From this, universality 
can only be deduced if all energy scales are sent to zero 
(see above). A fixed k — kp 0 sets a scale which be¬ 
comes active for all generic V{q) with V^^\q = 0) 7 ^ 0 
for some n G N, where V^^\q) denotes the n-th deriva¬ 
tive. Thus p^{k,uj) cannot expected to be universal on 
general grounds. In fact, a finite scale will destroy the 
scale invariance of the model and thus ‘quantum critical’ 
power-law scaling. In the case of the at g = 0 ‘infinitely 
fiat’ box potential this mechanism is not active as long 
as \k — kp\ < gc- 


We here supplement Ref. |9] and provide very strong 
evidence that the spectral function of the TLM at fixed 
k — kp 7 ^ 0 and for a generic potential is not char¬ 
acterized by power laws. The latter are only found if 
k — kp^O. To guide the reader we should from the out¬ 
set be very precise about our understanding of ‘power-law 
scaling’. We say that some real function / shows power- 
law scaling with exponent ^ close to xq G M from above if 
din \f{x)\/dln{x — xq) approaches ^ for x ^ xj. A sim¬ 
ilar definition can be given for power-law scaling from 
below. This does of course not exclude that / can to 
some degree be approximated by a power law or ‘resem¬ 
bles’ a power law for some range of x (close to xq) even 
if it does not fulfill the above criterion. We will return 
to this in Sects. liTTEl and |HIF| in which we present our 
results for p^{k,uj). 


More generally, we show that even for \k — kp\ ^ gc 
the spectral line shape depends on the details of V{q). 


Although the present study might be regarded as some¬ 
what technical - or even pedantic given that issues of the 
momentum dependence of V (g) in the TLM are virtually 
always nonchalantly ignored - our results have far reach¬ 
ing consequences. 
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C. Implications of our results 

The first implication of our results concerns the con¬ 
cept of LL universality. While p^{oj) of any model from 
the LL universality class shows the power-law suppres¬ 
sion of spectral weight for cj ^ 0 and p^(/cF, cj) a power- 
law peak or suppression at the same energy, LL univer¬ 
sality does not prediet power laws in p^{k^uj) for fixed 
k — ^ t). Evidently, if this type of universality cannot 

be proven in the low-energy fixed point model, the TLM, 
it cannot be a characteristic feature of the LL universality 
class. This does of course not exclude that eertain mod¬ 
els from the LL universality class might show such power 
laws, however, if so for more speeifie reasons than LL uni¬ 
versality. An obvious example for this is the TLM with 
box potential.!^ Other examples might be based upon the 
restriction of the (equilibrium) dynamics encountered in 
specific Id models with an extensive number of local con¬ 
served quantities (e.g. the Hubbard model) which are 
often Bet he ansatz solvable .1^^ A detailed discussion in 
which we relate our results to spectra of Id lattice mod¬ 
els obtained by either analytical or numerical approaches 
is given in Sect. |IV A[ 

Strongly linked to this are the implications of our find¬ 
ings for recent results on p^{k, uj) taking the nonlinearity 
of the single-particle dispersion in to accounfP^MUl which 
are embedded in the framework of the so-called ‘non¬ 
linear LL phenomenology ’.1^ In this power laws are not 
viewed as originating from quantum critical scale invari¬ 
ance but rather as resulting from a Fermi edge singu¬ 
larity like effect. In the phenomenological construction 
of the effective field theory including curvature effects of 
the dispersion the above described ad hoc regularization 
is employed. The spectral function is computed for this 
field theory. Our results obtained for linear LL theory 
raise the question whether the power laws found in ‘non¬ 
linear LL phenomenology’ are robust against a curvature 
of the bulk two-particle potential. We emphasize that 
these power laws are specific to the nonlinearity of the 
dispersion, which e.g. leads to momentum dependent ex¬ 
ponents, and are thus different from the nonanalyticities 
found for the TLM treated within the ad hoc procedure 
(or, for that matter, the TLM with box potential). Again 
this does of course not exclude that for speeifie Id mod¬ 
els, e.g. Bethe ansatz solvable lattice models, power laws 
with momentum dependent exponents might be realized. 
More on this is can be found in Sect. HYU 

Finally, our findings are of importance for the inter¬ 
pretation of experimental momentum resolved spectra. 
Even after decades of research none of the photoemis¬ 
sion experiments on quasi Id metals reporting on the 
observation of dispersing spin and charge peaks remains 
unchallenged.One reason for this is that, when inter¬ 
preting experimental data in the light of LL physics, the 
momentum resolved spectral function obtained within 
the TLM employing an ad hoc regularization is taken 
paradigmatically. Crucially, we find spin and charge 
peaks for generic V{q) and a sufficiently small ampli¬ 


tude of the two-particle interaction even though they are 
not given by power-law singularities. Our results show, 
however, that details are model dependent (in our case 
V{q) dependent) and therefore nonuniversal. Thus the 
detailed spectral features of the ad hoc regularized TLM 
cannot expected to be found in experimental spectra. 
Further account of the relevance of our results for ex¬ 


perimental spectra is given in Sect. IV C 


D. Structure of the paper 


The rest of our paper is structured as follows. In 
Sect. [ni we introduce the TLM and its bosonization solu¬ 
tion. Constructive bosonization of the field operator can 
be used to derive a closed analytical expression for the 
single-particle Green function, which, however, contains 
a momentum integral on the right hand side. This is dis¬ 
cussed in Sect. rmw In Sect, irir^ we introduce different 
versions of the ad hoc regularization to obtain int egral 


me 


and 


HID 


free expressions for G<{x,t). Sections 
are devoted to the technical details of how we obtain ex¬ 
act spectra for a box potential and arbitrary potentials, 
respectively. In Sect. ^we present results for the spec¬ 
tral function of the so-called 5 f 4 -model with intra-branch 
interaction only and different shapes of momentum de¬ 
pendency of the interaction. These are compared to the 
ones obtained by the ad hoc procedure. Section [III F| ist 
devoted to the spectral function of the spinless TLM - 
the spinless g 2 -gA-ioaode\. In Sect. |IV| we discuss the impli¬ 
cations of our results. When alluding to the relevance of 
our insights for the interpretation of photoemission data 
we in addition present spectral functions for the spinful 
TLM. 


II. THE TOMONAGA-LUTTINGER MODEL 


We here do not introduce the TLM by ‘deriving’ it 
from the Id interacting electron gas under certain as¬ 
sumptions (e.g. on the real-space range of the inter¬ 
action; no 2 /cf two-particle scattering processes^P or as 
the effective field theory for microscopic lattice modelJ^ 
but rather take it as a stand-alone model. It consists 
of independent right- {a = +) and left-moving (a = —) 
fermions with spin s, creation operators ctl. ^ disper¬ 
sion ^a{k) = av^ik — akp), density operators {q ^ 0) 
P<x,s{q) = T.kd,a,sH+q,a.,s^ ^nd particle number opera¬ 
tors nk,a,s = Following LuttingeP an infi¬ 

nite ‘Dirac sea’ filled in the ground state is assumed and 
thus the momentum quantum number k of both par¬ 
ticle species is unbounded. This simplifies the mathe¬ 
matical treatment as certain relations become operator 
identitieJ^ and are not only restricted to the low-energy 
part of the Hilbert space as in Tomonagas approach.^ 
This addition of states often requires normal ordering. 
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The Hamiltonian for a system of length L is given by 


H — ^ ^ ^ ^ ^OL (^) \j^k,OL,s {'^k,a,s) q\ 

k a,s 

+ ^ [9i,\\(.<i)^s,8' + 94,±{<i)^B-8'] PaA<i)pi,s'(.9) 

q^O ^ 

+ [92,\\{q)S8,s' + 92,±{9 )^s-s'] p+,s{9)P-,s'i9)-i'^) 

q^O 

s,s' 


Here (•. .)o denotes the (noninteracting) ground state ex¬ 
pectation value (normal ordering). We keep the explicit 
g'-dependence of the two-particle potential. The inter- 
{g 2 ) and intra-branch (^ 4 ) potentials are not necessarily 
equal and replace the potential V{q) referred to in the 
Introduction. Similarly, the interaction of spin parallel 
(II) and anti-parallel (T; this is a confusing but standard 
notatioiP) particles is not necessarily the same. If the 
TLM is considered as the low-energy fixed point model of 
the LL universality class this flexibility is required. The 
low-energy physics of any model from this class is charac¬ 
terized by four independent numbers, e.g. the two LL pa¬ 
rameters Kc/s and the two velocities Uc/s!^In the TLM 
for a given up (and k-p) those are fixed by the (in gen¬ 
eral) four independent ‘coupling constants’ = 0 ) 

(i = 2,4; hz =11,-L) at vanishing momentum transfer (see 
below). Therefore, to encounter nontrivial interaction ef¬ 
fects in the LL sense the = 0) should not all be 0. 

We restrict ourselves to these kind of interactions. 

We assume that the Fourier transforms ^i,«(g) of the 
two-particle interaction are even and for g' > 0 monotonic 
functions which vanish for \q\ ^ qc^ with an interaction 
cutoff gc- These requirements are physically sensible if 
the TLM is considered in its own right.^^If, however, the 
TLM is studied as the effective low energy model it is less 
clear if this assumption holds. It is thus crucial that it 
the assumption is not essential for our main conclusions. 
Relaxing it would merely complicate the calculations as 
positive and negative momenta would have to be treated 
seperately. We emphasize that at no stage of the discus¬ 
sion it will be necessary to introduce any further ultra¬ 
violet cutoffs ‘by hand’ despite the infinite (filled) Dirac 
sea at negative energies. In this sense the Hamiltonian 
Eq. 0 represents a mathematically well defined model. 

To be more precise the Hamiltonian in Eq. 0 defines 
a whole class of models as the four coupling functions 
di.nio) can be arbitrarily chosen as long as the introduced 
requirements are fulfilled. Still, we continue to refer to 
this class as the TLM. Note that a Hamiltonian of the 
form Eq. Q but with coupling constants gi^^ instead of 
coupling functions gi^^iQ) can often be found in the lit¬ 
erature. The necessary ultraviolet regularization is then 
left implicit and frequently not uniquely defined. 

We note that particle number contributions to the 
Hamiltonian which appear if the model is derived from 
the interacting Id electron gaJ^ are suppressed as they 
do not play any role for our considerations. 


The spinless version of the TLM follows from Eq. 0 
by dropping the spin index and keeping only gi (g) instead 
of for ^ = 2,4. 

Bosonization of the Hamiltonian and a canonical trans¬ 
formation lead tcPE 

^ = EE ‘^u{q)Pl{q)(^A9) + Eo, (2) 

q^O ip=c,s 


with bosonic operators ldA{q) describing collective 
charge (z/ = c) and spin {n = s) excitations (spin-charge 
separation) as well as the ground state energy Eq. The 
jSA(q) are linearly related to the densities pa,s{Q) of the 
fermions. The mode energies ujjy{q) are given by 


^u{q) 




Vu{q)A) 


where we have introduced the renormalized momentum- 
dependent charge and spin density velocities Vjy{q) and 
interactions 


9i,c/s{q) = [9i,\\{q) ±9i,±{q)] / 2 - ( 4 ) 

We already now emphasize that for momentum depen¬ 
dent giAo) f^o eigenmode dispersions will become non¬ 
linear. The linearization of the latter, that is the replace¬ 
ment 


Vv{q) v^{0) = v^, (5) 

is the crucial step in the ad hoc regularization procedure 
to derive integral-free expressions for correlation func¬ 
tions (see below). 

In Eq. 0 we have introduced the g-independent 
charge and spin velocities Vc/s relevant for the low-energy 
physics (all energy scales sent to 0) in the LL sense. The 
corresponding LL parameters K^/s of the TLM are ob¬ 
tained from 


KM = 


'1 +g4,t/(g)/(7ri>F) - 92,,y{q)/MF) 

1 + 9i,u{q)l Mf) + 92,v{q)l ( ttvf ) 


(6) 


in the limit g 0. Note that Kj^{q) = 1 if the inter¬ 
branch interaction g 2 ,v{q) vanishes. 

The spinless version of the bosonized Hamiltonian is 
obtained after dropping the terms with index v = s. 

When presenting our results for p^{k,uj) we will ini¬ 
tially focus on two special cases. The first is the spinful 
TLM with intra-branch interaction only, that is g2A^) ~ 
0, which is commonly referred to as the g 4 -model. The 
second is the spinless g 2 -^ 4 -model. They are paradig¬ 
matic for the two interaction effects characteristic for 
LLs: spin-charge separation and power-law scaling with 
interaction dependent exponents, respectively. Proceed¬ 
ing this way increases the transparency of our analysis. 
The scenario for the general spinful model can be deduced 
by combining the insights of both cases; in Sect. |IV C| we 
in addition present a few explicit results for the spinful 
TLM. 
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III. THE MOMENTUM RESOLVED SPECTRAL 
FUNCTION 

A. The single-particle Green function 

We are interested in the momentum resolved spectral 
function. It can be computed by Fourier transforms from 
the greater and lesser single-particle Green function 


The interaction enters in 


-j _ oo -| 

F{x,t) =2 X] H 


u=c,s n=l 


+ 2 'yt,{qn) (cos{qnx)e ‘“-('J’*)* - , (13) 


with 


= (V’a,5(a^,i)^/’i,s(0,0)) (7) 

iGa,s{x,t) = (V’t5(0,0)i/'a,s(a:,t)). (8) 

The field operators are related in the 

usual way to the creation and annihilation operators in 
momentum space 

■*/’«,«(a;) = -b X] (9) 

k 

Particle-hole symmetry of the TLM ensureJ^^^ 


= (14) 

and Qn = n2'KjL (periodic boundary conditions). Due to 
the decay of the on the scale the momentum 

sum in F(x, t) is convergent in the ultraviolet. Note in 
particular that for ^ qc the two terms in the first line 
cancel each other as ujy{qn) '^F^n in this limit [com¬ 
pare Eq. ©]• The term in the second line of Eq. is 
convergent as ^y{qn) 0 for q^ ^ qc [compare Eqs. (|14|) 
and (p)]. Eor vanishing interaction F{x^t) = 0. Thus 
[G^] (x, t) is the noninteracting Green function. In the 
thermodynamic limit L ^ oo it becomes 


Ga,s{x,t) = G<^,{-x,-t). (10) 

and it is sufficient to consider the greater Green function. 
We note in passing that band filling is not an issue in the 
TLM as we consider it. However, when the model is in¬ 
vestigated as the effective low-energy model of another 
Id correlated electron model the band filling of the latter 
will enter in the LL parameters and velocities character¬ 
izing the low-energy physics.l^ 

Eurthermore, with x —x we can go over from the 
Green function of right-movers to the one of left-movers. 
Therefore, we will only study G^ s(^ 5 ^)- fhe absence 
of a magnetic field the Green function is spin independent 
and we thus suppress the spin index from now on. 

To compute the ground state expectation value in 
Eq. 0 we use constructive bosonization of the field 
operator .1^ We emphasize that this approach is based 
on operator identities and does not require the introduc¬ 
tion of any cutoffs if one first considers finite systems 
of length L (as we will do). In more phenomenological 
approacheJ^a cutoff is introduced - often denoted hy 1/a 
and referred to as an ‘effective band width’ - which for¬ 
mally has to be send to infinity. However, it is frequently 
kept finite artificially. This is part of one of the possible 
ad hoc ultraviolet regularizations. 

The exact greater Green function of right-movers for 
the most general spinful g 2 -gA-m.ode\ is given by 

iG^{x,t)=i (11) 

with 

i{kF+F/L)x 

[G^] {x,t) = - - - -exp -F‘i-G-VEt+^o+) I ^ 

( 12 ) 


[GV\yx,t)=i 




27r X — v^t -h i0+ ’ 


(15) 


The factor exp (—gnO^) in Eq. (12) ensuring convergence 
appears naturally, and is not related to any ad hoc reg¬ 
ularization. Onl y w ith this factor the exactly known 
[G^] (x, t) Eq. ilbl and from this the exact noninter¬ 
acting spectral function 


[pX{k,u})]° {x,t) = e{k - kF)S{u} - ^+{k)) (16) 

can be obtained. Due to the linear single-particle dis¬ 
persion [G^] (x, t) is of relativistic form. The greater 
spectral function is defined as 

1 

p^{kn,uj) = — / / dx iG^{x^t). 

2tt J-L/2 

_ _ (17) 

As G^(x,t) Eqs. (Im to is an analytic function in 
the lower half of the complex t-plane p^{k^uj) has non¬ 
vanishing weight only for cj > 0. 

To compute p^{k,uj) for arbitrary potentials in the 
thermodynamic limit three nested integrals have to be 
performed. The integrands are slowly decaying, oscilla¬ 
tory, and have poles close to the real axis. Therefore, the 
accuracy which can be achieved when straightforwardly 
performing the integrals numerically, given gi,/^(g), is not 
sufficient to answer the question whether or not power 
laws can be found for k — k-p ^ 0. We are thus forced 
to proceed differently. Before presenting our approach 
to the exact spectral function of the TLM in Sects. [IH^ 
and |HID] we will describe approximate ad hoc procedures 
which were pursued in the literature to make analytical 
progress. 
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B. The ad hoc regularization 

In the natural way of writing the exact Green function 
of the TLM Eq. 0 the noninteracting one is factor¬ 
ized out P However, othe r express ions for can 

be found in the literatureThese result from the 
following procedure. The term in the curly brackets of 
Eq. (12) is canceled against the last term in the first 
line of Eq. (13). After this step the remaining q-sum in 


Eq. (13) is no longer convergent in the ultraviolet. Two 


ways have been reported how to deal with this problem. 

In the first (i) the remaining first term of the first 
line of Eq. (13) is multiplied by exp(—To ob¬ 
tain integral-free expressions for one then as¬ 

sumes that 7zy(gn) is given by exjp {—QnA) which can 
be reached by choosing a proper momentum dependence 
of the gi,K{q). Thus one simply selects a certain interac¬ 
tion potential. In addition, one linearizes the eigenmode 
dispersions uJu{q) v^q for all q. This is not a conse¬ 
quence of the special choice of the interaction potential 
but is rather an approximation done independently ‘by 
hand’. In the thermodynamic limit one then obtains 


oik-px 


iGiux.t)=— n 


1 


nl/2 


X — Vyt -h i0 + 


[x — Vjyt + iA) {x + Vyt — iA) _ 


7G2 


(18) 


In the second ad hoc procedure (ii) the entire (remain¬ 
ing) argument of the sum in Eq. (13) is multiplied by 
exp (—g^A) with a finite momentum cutoff 1/A > 0. The 
momentum dependence of "fu{Qn) is dropped and uJu{q) 
is again linearized. This leads to 


f^ikpx 


n 


1 


nl/2 


X — Vyt + iA 
A^ 


(x — Vyt + iA) [x -h Vyt — iA) 


7G2 


(19) 


Einally, a third way (iii) to obtain an integral-free ex¬ 
pression for G^(x,t) is based directly on Eqs. (11) to 
In this the momentum dependence of 7zy(^) is 
dropped and ujy{q) is linearized after the square bracket 
in Eq. (13) was multiplied by exp (—g^A). This way one 
obtains 


[G^](iii)(x,t) =- 


pikpx 


1 


2tt x — v^t + i0+ 


n 


v=c,s 


x — v^t + iA 
x — Vyt + iA 

A^ 


1/2 


17^2 


{x — Vyt -h iA) {x -h Vyt — iA) 


(20) 


We note that for the special case 7 s = 0 (spin-rotational 
invariant interaction) and Vs = v-p this is exactly the 
expression derived by Dzyaloshinskii and LarkirP^ in a 
purely fermionic approach which is based on Ward iden¬ 
tities and the closed loop theorem.l^ 

Obviously, all three approximate functions differ. How¬ 
ever, their asymptotic behavior (power-law decay) for 
large space-time arguments is the same, in particular 
along the four special directions x = ±.Vyt. Eor this rea¬ 
son they all lead to the same power laws in pf_{k,uj) for 
uj close to ^Vy{k — k^) as derived in Refs. [5l [9l [131 

M 

Similarly, for x = 0 all three approximate expres¬ 
sions have the same asymptotic behavior at large |t|, 
^ with = 7 c + 7 s (and are analytic functions 

in the lower half of the complex t-plane). After a sin¬ 
gle Eourier transform from time to frequency this gives 
p^(cj) ^ Q{uj)uj^. As discussed in Ref. [9] this power-law 
scaling of the momentum integrated spectral function is 
also found based on the exact expressions (11) to (13) 
and for arbitrary gi^^ig) (fulfilling the above mentioned 
restrictions) with jy ly{fL — 0 )- This universality is 
based on the RG irrelevance of the momentum depen¬ 
dence of the interactional the appearance of the power 
law with an exponent set by the two-particle potential 
at vanishing momentum transfer g = 0 is not affected by 
the potential away from this point. Any of the discussed 
ad hoc procedures (i) to (iii) can thus be employed with¬ 
out spoiling the universal behavior. The same holds at 
t = 0 but X 7 ^ 0 which after Eourier transform leads to 
the momentum distribution function n+(/c) which also 
shows universal power-law scaling for k The 

question we address here is whether or not the ad hoc 
regularized Green functions can also be used to obtain 
universal results for p^{k,uj) at k — k^ 0. 

As already emphasized the linearization of the ujy{q) 
is the crucial step which leads to power laws in p^{k,uj) 
for k — k^ 7 ^ 0 in the ad hoc procedures. In Ref. [9] it 
was shown using a generalization of the stationary phase 
method that any g = 0-curvature of Vy{q) affects the 
asymptotic behavior of G^(x,t) in the distinguished di¬ 
rections X = :Avyt. The decay is no longer solely given by 
the ^y and therefore not only by the ^i,/^(0). However, in 
the ad hoc procedures (with constant Vy), one can ana¬ 
lytically show that it is the asymptotic behavior in these 
directions of the x-t-plane which leads to the power laws 
at UJ = -^Vy[k — /cp]. This raises doubts that the latter 
are generic. 

Besides the spinful TLM we also consider its spinless 
version. Erom the above expressions for G^(x, t) [includ¬ 
ing the ones of the ad hoc regularization (i) to (iii)] the 
spinless Green function is obtained by setting 7 ^ = 7 ^, 
ujc = uJs, and Vc = Vg. All what was said in the last 
two paragraphs about power-law behavior and universal¬ 
ity remains valid in the spinless case up to the (obvious) 
exception that in the ad hoc regularized spectral func¬ 
tion only the two (instead of four) distinguished energies 
:Avc{k — k-p) exist. 
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We note in passing that it was realized decades ago 
that within the ad hoc procedures (i) and (ii) exact spec¬ 
tral sum rules are not fulfilled 


C. How to compute p^(k^uj) for a box potential 


In the Introduction we mentioned that exact results 
for the momentum resolved spectral function of the TLM 
were derived based on Eqs. to ( [l^ assuming a box 
potential!^ We will compare our results for other po¬ 
tentials to these. To be self-contained we here give all 
the formulas required to obtain p^{k^uj) of the spin¬ 
ful ^ 4 -model and the spinless ^ 2 -^ 4 -mc)del for gi^,^{q) = 
9i,K^{qc ~ More details are presented for the first 
case which was not separately discussed in Ref. [151 We 
note that for a reader primarily interested in results it is 
not necessary to understand the technical details of this 
section in full detail. 


1. The spinful gA-model 


For 5 ^ 2 ,/t(^) = 0 it directly follows that 7zy(g) = 0 as 
Kjy(^ = 1; compare Eqs. (14) and •HI- Equations ( pT] ) 
to (13) for the exact Green function thus simplify consid¬ 


erably. The same holds for the ad hoc regularized expres¬ 
sions Eqs. (18) to (20) as = 0. These are characterized 
by square-root singularities at x = Vct and x = Vgt. After 
Fourier transform they lead to square root singul arities 
in p^{k,uj) for uj ^ Vc[k — k^] and uj ^ Vs[k — 

We note that [G^](iii)(x, t) contains the additional factor 
{x — v-pt + iA) l{x — v^t + i0+) which does not drop out 
as A > 0. Further down we will discuss how this term 
affects the spectral properties. As the ground state of 
the 5 f 4 -model remains the noninteracting onef^ p^{k^uj) 
has finite weight only for /c > /cp. 

For the case of a box potential, it is 


Vv{q) = 


Vu, q<qc 

VF, q > qc- 


( 21 ) 


Setting for convenience 2 : = exp{i(27r/I/)(x — v-pt)} as 
well as Zjy = exp{i(27r/I/)(x — Vj^t)} and expanding the 


exponential function one obtains from Eqs. 0 to @ 


iG>{x,t) = le*(2-/i)(«F+i)^ I g , 


kI=0 




n=l \m=0 
/ 00 


m\ 


(l/(2n))-^ ^ 


n] 


0=0 


j! 


E 

v/=o 


(l/(2n))y„, 


( 22 ) 


_L i(27r/L)(nF+l)a: 

L 


E 

Vm=0 





(''^c) I \ 

i 5 


(23) 


with Qc = ncfl'KjL. Here k^ = (2nF + l)7r/T and np is 
the index of the last occupied single-particle state; k^ 
thus lies in between the last occupied and the first un¬ 
occupied one. The coefficients in Eq. (23) can be deter¬ 
mined by a recursion relation, where m > 1, / G No and 
i = 0,..., m — 1: 


Am) 


= E 

j=0 


( XjTfl)^ (m—1) 


Am) _ ST' (l/(2m))-^ Am-i) 
^Im+i / j p m(Z— 

j=o 


= E 

j=0 

The initial values are given by 


.(1) 


= E 


i-iy 


(24) 

(25) 

(26) 


j=0 


and = (1/2)^/ml. The recursion can easily be per¬ 
formed on a computer. The double Fourier transform 
can be computed analytically and one obtains 


n n—l 

1=0 i=0 

2tt 

X 5 uj - - [{n-l - j)vF + Ivc + jVs) 


(27) 


with n = n — (np + 1). This way the exact spectral func¬ 
tion p^{kn,uj) of the TLM with box potential ca n easil y 
be computed for large but finite systems (see Sect. HIE). 

To obtain analytical insights we rewrite the Green 
function as 

_ ^ i(27r/L)(nF + l)x 


iG^{x, t) = — e 


yin 


V n=l ^ ^ n=nc+l 


(28) 





















Thus 


2 . The spinless g2-g4-model 


E 

m=0 


CX3 

1+ E 

n=nc + l 


n 


1 

2 


E 


Vn=nc + 1 



+ ... 


(29) 

and we immediately see that = 1, = 0 for 

1 < m < Tie and = 1/m for ric + 1 < m < 2nc + 1. 
For n < Tie the simplified expression 


CJ - ^ {hVs + l{Vc - Vs)) 

(30) 

holds. For fixed ft and Vg < Uc, there is only spectral 
weight for Vs[k — k^] < (v < Vc[k — k^] (up to corrections 
of order 1/L). For Vg > Vc the roles of the two velocities 
are interchanged. We note that for k — k-p < Qc the bare 
Fermi velocity up drops out. 

Further analytical results can be obtained employing 



where ( ^^) is the generalized binomial coefficient. From 
this it follows that = 1 and 

^ const. X j- 1/2^ (32) 


where 

(/) ^ ^ ^ 

was used. Inserting this for energies close to Vg[kn — k-p] 
into Eq. ( [M| ) , where I is the integer next to {uj — Vg [kn — 
k^]) /[{27rjL){vc — Us)], we obtain for L ^ oo the one¬ 
sided square-root singularity {vc > Vg) 

p^{k,uj) ^ 0 (cj — Vg[k — /cp]) (<^ — Vg[k — /cp])~^^^ . 

(34) 

Analogously we get close to Vc[kn — k-p] 

p^{k,uj) ^ 0 (cj — Vc[k — /cp]) {uj — Vc[k — /cp])~^^^ . 

(35) 

For Vg > Vc we only need to interchange the two veloci¬ 
ties. 

We have thus shown that the exact spectral function 
of the ^ 4 -model with box potential shows the edge sin¬ 
gularities also found within the ad hoc procedures (i) to 

(iii). 

We note in passing that for k = k-p the spectral func¬ 
tion reduces to a (^-function of weight 1 located at cj = 0. 


pXikn.uj) = }_^bi 


hMAric) 


'^h—l 


In Sect. IV of Ref. [15] a recursive way of computing 
the spectral function for the full spinless TLM with box 
potential similar to Eqs. (24)-(27) was introduced. It is 
given by 


n+r 


r=max{0,—n} 1=0 

27r 

X 6 uj -— ((h + r — 1)vy + (r -1- l)vc) 


(36) 


with A = exp{J]]^^^(l/n)}. The coefficients are de¬ 
termined as in Eq. and ( [^ . Eor the recursion 
relation reads 


Am) f 

im+i j! I rn 

3=0 ^ 




(37) 


with the initial values 6^^ 


„(”c) 


= (1 -h jc)'^l'm\, and for the 


Jni) (m-1) 

4m+i -I m(Z-i)+i’ 

3=0 


(38) 


with the initial values Cm^ = jml. 

In analogy to the considerations for the spinful g^- 
model, we can infer from the behavior of the coefficients 
the behavior of the spectral function close to Vc[k — k^] 
(for L oo). Eor k — k^ > 0 we fineP^ 

p^{k,u!) Q{u!-Vc[k-kF]) {u - Vc[k - kF]y°~'^, (39) 

that is a divergence if 7 c < 1 (not to strong interactions). 
In contrast to the ^ 4 -model p^ can now also carry spec¬ 
tral weight for (k — k-p) < 0. At the threshold we obtairP^ 

pX{k,uj) ^ 0(-u; - Vc[k - /cp]) {-u - Vc[k - /cp])^" , 

(40) 

that is a power-law suppression since 7 c > 0 in the full 
model. In the special case k = k-p, power-law behav¬ 
ior with p^{k-F,uj) oc is found. These threshold 

power laws can also be found based on the Green func¬ 
tion of the ad hoc procedures (i) to (iii) discussed in 

Sect.lmlpniEHni 


D. How to compute p^(k,(jj) for arbitrary potentials 

We now show that expressions for the Green and the 
spectral function which involve recursively computed co¬ 
efficients can also be given for an arbitrary momentum 
dependence of the two-particle potential. Eor the box po¬ 
tential the dispersion of the elementary charge and spin 
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modes is piecewise linear. This changes for arbitrary po¬ 
tentials. As a consequence the coefficients become time- 
dependent and the Fourier transform with respect to time 
has to be performed numerically, e.g. using fast Fourier 
transform (FFT). 


1. The spinful g^-model 


For arbitrary potentials the spin and charge velocity 
are no longer piecewise momentum independent. If the 
potential is effectively zero for q > with a Qc which 
we take sufficiently larger then the characteristic scale 
Qc^ we can work with g 4 ,^y{q) 94,y{Q) ^{Qc ~ q‘^) 

all practical purposes. Then we can rewrite T(x, t) in 
Eq. p) 


The remaining Fourier transform 

/ oo 

dt iG\{kn,t) (47) 

-OO 

can be performed numerically as a FFT. Since for the 
finite system the spectral function consists of (5-peaks, 
the Green function does not decay in time. Therefore, 
we have to multiply with a decaying func¬ 

tion before performing the FFT. Here, we will always 
use the exponential function exp{—x|t|}. In frequency 
space, this corresponds to a convolution of the spectral 
function with the Lorentzian 'K~^xj{uP‘ + y^), i.e. each 
(5-peak is broadened into a Lorentzian of width y. 


2. The spinless g 2 -g 4 -Tnodel 




n=l \ 

g-i(27r/L)nvFt^ ^n 


(41) 


with 2 ) = exp{i(27r/I/)x}, and use time-dependent coeffi¬ 
cients to write 

iGX{x,t) = i e*(Wi)("F+l)a; 


ric / oo 

HE 

n=l \ m=0 


(l/») 

ml 


1 


— iuJc(qn)t 


_ g-i(27r/L)nvFt 


(42) 


^ _ gi( 2 ,r/L)(np + l):r ^ ^(n,) (43^ 


L 


m=0 


The recursion relation for the (t) is given by 




^ I” —i{27T I L)vFt ^ 


1=0 


\^-iuJc{qi)t \^-iuJs{qi)t _ ^-i{2T^jL)vFt 


n I 


(44) 


a 


(m) 


<«) = E 


(1/m)-^ 


j=0 


\^-iuJc{qm)t 


^ f^-iuJs{qm)t _ g-i(27r/L)mvFt 


a 


(m-l) 


(45) 


In the same way as for the spinful ^ 4 -model, we can 
introduce recursively defined time-dependent coefficients 
for the spinless ^ 2 -^ 4 -inodel 


a, 


" ^ 1 — 7 

d)(t) = ^ - p-h27r/L)7;Ft] 


1=0 


X [(1 7c(?i)) 




j=0 

_ ^—i{2T^lL)ravFt 




a 


(m-l) 


it) 


and 




(48) 


(49) 


(50) 


= E 

j=0 


{1/my 


iMmY 


—iujc(qm)t 


^ (t) 


m{l-j)+i 


(51) 


2 fc{Qn) 


n 


In terms of these we can rewrite 

1 r 

iG^{kn,t) = -2^] 

f n=l j 

00 

X E ^n+liObi^^Ht) (52) 

r=max{0,—n} 

and the remaining Fourier transform to obtain p^{kn,oj) 
can be performed, again after multiplication with 
exp{—y|t|}, numerically by means of a FFT. 


where in the second line m > 1, / G No and i = 0,..., m — 
1. With this representation of the Green function, the 
Fourier transform to momentum space can be performed 
analytically 

iG^{kn,t) = (46) 


E. Spectra of the spinful ^4-model 


Based on the formulas given in Sects. [HI C l| and pll D 1 


we are in a position to compute the exact spectral func¬ 
tion of the spinful ^ 4 -model for different potentials at 
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finite system size L. For the box potential both Fourier 
transforms can be performed analytically while for arbi¬ 
trary potentials the time-transform is performed numer¬ 
ically as a FFT. In this case only broadened spectra can 
be obtained. We will compare the results to those derived 
from one of the ad hoc procedures of Sect. Enl 
Besides the box potential we consider {v = c, s) 


54^(5)= 54,. exp{-(g/9c)V9} 

(P=4), 

(53) 

5fr(</)=54,.exp{-(<//9c)"} 

(gauss), 

(54) 

= 54,.exp{-3 \q/qc\} 

(exp). 

(55) 


The factors in the exponential function were chosen such 
that, besides at g = 0, all potentials have the same value 
at qc = 3gc; at this momentum they have decayed to 
~ 10“^ of the q = 0 value and we can safely set the 
potentials to 0 for g > gc- Considering larger qc we have 
verified that this does indeed not affect our results. For 
small momenta the potentials go as 1 — gA,v{(l)/Q a.v ^ 
Ig/gd^, with p = oo for the box, p = 4 for the ‘p=4’- 
potential, p = 2 for the Gaussian potential, and p = 1 for 
the exponential potential. The exponent p is a measure 
for the flatness of the potential at g = 0 ; see also Ref. [9l 

For definiteness in our calculations we have always cho¬ 
sen p4,c/('7rt’F) = 1/2 = —gA,s/{'T:vY) without affecting 
our conclusions. In this case Vc > Vg. The measure for 
the system size is ric- To numerically compute the re¬ 
cursively defined coefficients within reasonable time we 
choose Uc = 5-10^. The broadening y is chosen such 
that in the broadened spectral function no effects of the 
individual J-peaks are visible. For the given Uc we take 

x/{vFqc) =5 • log'. 

In Figs.l^) to^) we show p^{k,uj) for three different 
small {k —The curves labeled as ‘box (ana)’ are 
the results for the box potential obtained from Eq. ( [M| ) 
without convoluting it with a Lorentzian. The weights 
of the (5-peaks are divided by the level spacing (which 
for the given parameters is l/ric) and are connected to 
form a smooth curve. All other results (up to ‘broadened 
(iii)’; see below) were obtained from Eq. ( [4^ multiplied 
with an exponentially decaying function exp{—x|t|} and 
transformed with a EET; for comparison we also show 
the broadened spectra for the box potential [box (EET)]. 

As can be seen from Eig. Ef), for k very close to k-p 
all curves are nearly indistinguishable. The broaden¬ 
ing of the spectra is visible by the weight ‘leaking out’ 
for uj < Vs[k — k^] and uj > Vc[k — /cf]. However, al¬ 
ready for {k — /cF)/gc = 0.1 [Eig. [^)], the spectral func¬ 
tion calculated with an exponentially decaying potential 
shows pronounced differences to the other ones. Eor 
{k — kp)/qc = 0.3, also the curve of the Gaussian po¬ 
tential deviates from the one of the box potential, see 
Fig. [^). The spectrum of the ‘very flat’ p=4 potential 
still lies on top of the one obtained for the box potential; 
differences only appear at even larger k — kp (not shown). 
Obviously, the smaller the above introduced index p, that 
is the less ‘flat’ the potential is at g = 0, the faster the 
line shape of the spectra starts to deviate from the one 





FIG. 1. (Color online) Spectral function of the spinful g 4 - 
model as a function of energy for a) (k — kY)/qc = 0.01, b) 
(k — /cF)/gc = 0.1, and c) (k — /cF)/gc = 0.3. Spectra for the 
different potentials are shown. In addition to the FFT data 
with broadening x/(r’Fgc) = 5 • 10“^ we show unbroadened 
results for the box potential [‘box (ana)’] and the broadened 
results of the ad hoc procedure (iii) [‘broadened (iii)’]. The 
other parameters are g4,c/(7ri’F) = 1/2 = —g4,s/(7ri’F) and 
Uc = 5 • 10^. The curves are partly hidden by others (see 
the text). The inset of c) shows the position of the spin peak 
extracted from the data (stars) in comparison to the collective 
spin mode dispersion ujs{k — k^) (lines). The dashed line 
displays the unrenormalized dispersion v^ik — /cf). 
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of the box potential when k — increases. 

The maxima of p^{k,uj) are located at ujjy{k — /cp), 
instead of at Vy[k — k-p]; for \k — k-^l < Qc both positions 
are equal for the box potential. This is shown in the 
inset of Fig. El), where the spin dispersion relation for 
the different potentials (full lines) is compared to the 
numerically determined maxima of the spectra (stars). 
The agreement is very good; the charge peak behaves 
similarly. In the limit of small k — k-p and large p the 
difference between uJu{k — kY) and Vulk — k^] is negligible. 

In the literature the third ad hoc regularization proce¬ 
dure leading to [G+](iii)(^, t) Eq. (20) is considered ‘the 
best’ on d^ l ^^ l as it leads to a [p^](iii)(/c, cj) which fulfills 
exact sum rules For = 0 a closed analyt¬ 

ical expression of the double Fourier transform can be 
given; see Eqs. (3.20) and (3.21) of Ref. [28l We therefore 
added a graph of this analytical result convoluted with a 
Lorentzian of width y as the dashed lines in Figs. Ef) 
to El) (with 1 /A —> Qc). As mentioned above these 
curves as well as the exact function for the box potential 
show single-sided (threshold) square-root singularities at 
v^lk — k-p] (for a more detailed analysis on this see below). 
The approximate spectrum [p^]qii)(/c, cj) is characterized 
by an additional feature at uj = v^lk — k^] (barely visible 
in Fig.[^) for {k — kY)/qc = 0.01); as analyzed in Ref. [28] 
a logarithmic divergence appears at this energy. It results 
from the additional factor {x — v-pt + iA) /{x — vpt + i 0 +) 
of [G^](iii)(x, t) Eq. (20). Obviously, the exact spectral 
function of the ^ 4 -mMel does not display this feature 
for any of the considered potentials. For the box poten¬ 
tial, which besides the missing peak at uj = v^lk — k-p] 
shows the same spectral characteristics as found in the 
ad hoc procedure, this can even be understood analyti¬ 
cally: as seen in Eq. ( [3Q| ) for k — k-p < Qc the information 
on vy drops completely out. The logarithmic divergence 
of [p^](iii)(A:,cj) at uj = VY[k — k^] is thus an artifact of 
the ad hoc regularization. As discussed very recently, this 
logarithmic divergence ioi k — kp ^ Qc turns into a power- 
law one if the spinful ( 02 -^ 4 ) TLM is treated within the 
ad hoc procedure (hi) .[^ Also this feature is an artifact of 
the ad hoc procedure. With Ref. [29] in mind we empha¬ 
size that this does not exclude that for /c — /cp ^ all the 
spectral weight is located around uj = VY{k — /cp), which, 
in fact, is generically the case. This was discussed for 
the box potential in Ref. [T5] and for general interactions 
in Ref. |27| We here do not investigate this any further 
as we are exclusively interested in the spectra at small 
\k-kY\. 

In order to investigate whether or not the spectra at 
fixed t) < k — kp qc and for uj close to the maxima 
at cjmax show powcr-law scaling we do not simply want 
to rely on the quality of power-law fits to the broadened 
data. Instead we take the logarithmic derivative 


logder(c(;) = 


din [p^{k,uj)] 
din \uj 

^max I 


(56) 


of our data with cjmax equal to the corresponding peak 
positions. It is numerically approximated by centered 


differences and provides a very sensitive approach in the 
search for power laws. If for uj cjmax, logder(cc;) ap¬ 
proaches a constant f the spectral function shows power- 
law scaling with exponent f close to cjmax (according to 
the definition given in the introduction). 




FIG. 2 . (Color online) Logarithmic derivative Eq. (56) of the 
spectral function of the spinful ^ 4 -niodel close to the spin peak 
for a) {k — kF)/qc = O.OI and b) (/c —/cf)/^c = O.I. In addition 
to the broadened FFT data for the different potentials and 
the unbroadened results for the box potential [‘box (ana)’] we 
show broadened [‘broadened pi’] as well as unbroadened [‘pi’] 
data for the product of single-sided square-root singularities 
Eq. (57). The parameters are as in Fig. 


In Figs.|2|i) [{k-kF)/qc = 0.01] and[^) [{k-kF)/qc = 
O.I] we show logder(u;) close to the spin peak for the dif¬ 
ferent two-particle potentials. For symmetry reasons the 
behavior close to the charge peak is the same. Instead of 
an ad hoc regularized spectral function - which is spoiled 
by the spurious peak at uj = VY{k — kp) - we this time 
present the logarithmic derivative of the simple normal¬ 
ized product of two single-sided square-root singularities 


[pX]pi{k,uj) 


I & {uj — Vs[k — kp]) & {vc[k — kp] — uj) 

^ {u!-Vs[k- kp])^^'^ {vc[k - kp] - ’ 

(57) 


for reference. It is indicated by ‘pi’. In addition we 
present the logarithmic derivative of this expression con- 
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voluted with a Lorentzian of width x, indicated by 
‘broadened ph. The unbroadened data for the box po¬ 
tential - for which we analytically know that for L ^ oo 
a square root singularity at Vslk — kp] exists [see Eq. (34)] 
- very nicely follow the ‘pi’ curve down to cj — cJmax ~ 
4- At this energy finite size eorreetions destroy 

the power-law scaling. This can be verified by considering 
different Uc, that is different system sizes. The ‘broad¬ 
ened pi’ curve starts to deviate from the unbroadened 
one at cj — c^max ~ 6 • From this we conclude 

that for the chosen and y the broadening y cuts off the 
power-law scaling at slightly larger energies than the sys¬ 
tem size. This is consistent with the observation that the 
broadened spectrum of the box potential [‘box (EFT)’] 
almost perfectly follows the ‘broadened pi’ curve down 
to the much smaller scale uj — c^max ~ 6 • IO^^'T’f^c- This 
gives US confidence that the EFT data for the other po¬ 
tentials are unspoiled by both finite size and broadening 
effects down to cj — c^max ~ 6 • IO^^'T’f^c- At the ‘high 
energy’ end possible power law scaling close to the spin 
peak is cutoff by the charge peak. The latter is the ori¬ 
gin of the dominant feature at cj — cjmax ~ and 

^ — ^max ~ 10 “^'r’F^C 5 respectively. 


While for {k — /cf)/^c = 0-01 the data of the ‘p=4’ 
and the ‘gauss’ potential lie perfectly on top of the ‘box 
(EFT)’ data, and one is tempted to conclude that they 
are consistent with a square-root singularity at cjmax, the 
data for the exponential potential clearly differ and are 
inconsistent with such behavior. For {k — /cf)/^c = 0-1 
in addition the data for the Gaussian potential are in¬ 
compatible with this type of singularity. For even larger 
k — k^ (but still smaller than ; not shown) also the data 
for the ‘p=4’ potential no longer follow the ones of the 
box potential. 


The most consistent interpretation of our results is that 
for any potential which is not ‘infinitely flat’ at g = 0, 
that is if p < oo, the spectral function at fixed k — k^ > 0 
does strictly speaking not show power-law scaling close to 
<^max- The less ‘flat’ the potential is, that is, the smaller 
p, the faster this becomes apparent as /c — /cf is increased. 
The square-root singularities found for the box potential 
(and for the ad hoc procedures) are cut off by the curva¬ 
ture of the potential close to g = 0. They can thus not be 
considered as universal features of the spinful p 4 -model. 

We note that the data for p < oc are not only incon¬ 
sistent with power-law scaling when taking cjmax as the 
point of reference. We studied the behavior relative to 
other distinguished energies (e.g. Vj^[k — k-p] and the - 
due to the broadening - apparent thresholds). For none 
of these we find behavior which is consistent with power 
laws. 


Despite this lack of power-law behavior, for all interac¬ 
tion potentials studied the exact spectral function of the 
g 4 -model is still characterized by spin and charge peaks. 


F. Spectra of the spinless g2-g4-niodel 


Employing the formulas given in Sects. Ill C 2 and 
|IIID 2| we can compute the exact spectral function of the 
spinless TLM - the spinless g 2 -g 4 -mode\ - for different 
forms of the potential. In the last section we saw that 
for the small k — k^ we are interested in, the spectra for 
the ‘p=4’ potential Eq. (53) barely differ from the ones 


obtained for the box potential. The same holds for the 
full spinless model and in this section we focus on the 
box potential, the Gaussian one Eq. (54) as well as the 


exponential potential Eq. (55). To prevent an inflation of 


cases and parameters we consider the physically reason¬ 
able situation of equal inter- and intra-branch scattering 
92 {q) = 94 ( 9 ) = 9 {q)- We choose g(0) such that Vc = 2 vy 

and 7c(0) = 1/8. 

To compute p^{k,uj) via G^(fe,t) Eq. ( [^ and EFT 
for general potentials or via Eq. (36) for the box po¬ 


tential we have to recursively compute more coefficients 
and perform additional sums in comparison to what was 
necessary in the spinless g 4 -model. This increases the 
numerical resources required and we thus have to con¬ 
sider smaller system sizes compared to the latter; we 
choose ric = 2 • 10^. Furthermore, at fixed Uc the en¬ 
ergy level spacing of the spinless g 2 - 9 A-iociode\ is larger 
than that of the spinful g 4 -model. To obtain smooth 
curves we thus have to increase the broadening y. We 
take y/('r’Fgc) = 10“^. 

In Figs.j^) to|^) we present the total speetral funetion 
p^{k,uj) = p^{k,uj) ^ p^{k,uj) for {k - kF)/qc = 0, 0.01, 
and 0.1. We switched to this as it simultaneously shows 
the photoemission as well as inverse photoemission part 
of the spectrum. For increasing k — k-p the photoemis¬ 
sion part looses weight quickly. Thus in Fig. [^) with 
{k — /cF)/gc = 0.1 we only present a zoom-in of the in¬ 
verse photoemission part. In addition to the broadened 
functions for the box, Gaussian, and exponential poten¬ 
tials obtained by EFT we show the unbroadened one of 
the box potential [see Eq. ([^]. As above, to obtain the 
latter the weights of the ^-peaks were divided by the level 
spacing (which for the given parameters is 4/nc) and are 
connected to form a smooth curve. 

As discussed in Sect. |III C^ for 0 < k — k^ < qc and the 
box potential p+(/c,Cc;) shows threshold power-law non- 
analyticities at dzVc[k — k-p] with exponents 7c — 1 (for 
cj > 0) and jc (for cj < 0). In ‘box (EFT)’ these are 
broadened. Similarly to the spinful g 4 -model we observe 
that the smaller p the faster the line shape starts to de¬ 
viate from the one of the box potential when increasing 
k — /cf- For {k — /cF)/gc = 0.1 [see Fig. m and the 
exponential potential with p = 1 this already leads to 
a strongly modified distribution of the spectral weight. 
The deformed line shape can be understood in due detail 
when coniparing it to the spectral function of the spinless 
^ 4 -model.l^ This detailed analysis is beyond the scope of 
the present paper. 

From Ref. [9] and the RG irrelevance of the momentum 
dependence of the two-particle potential we expect that 
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FIG. 3. (Color online) Total spectral function of the g 2 -g 4 - 
model as a function of energy for a) (/c — /cf)/^c = 0 , b) 
{k — kF)/qc = 0.01, and c) (/c — /cf)/^c = 0.1. Spectra for 
the different potentials are shown. In addition to the FFT 
data with broadening x/(r’F^c) = 10“^ we show unbroadened 
results for the box potential [‘box (ana)’]. The interaction at 
^ = 0 is chosen such that Vc = 2i’f and 7 c (0) = 1/8. The 
system size is set by ric = 2 • 10^. 


for all potentials p-^{k-F^uj) ^ On first glance 

the data of Fig. [^) appear to be consistent with this 


behavior, however, a more thorough analysis is required. 
In Fig. [^) we plot the logarithmic derivative Eq. (56) 
(with cjmax = 0) of the cj > 0 broadened FFT data of 
Fig#) as dotted lines. While the logarithmic derivative 
of the unbroadened ‘box (ana)’ data (solid line) nicely 
shows a plateau at the expected exponent 27c —1 = —3/4 
(for the given parameters), which is only spoiled at very 
small uj/{vFqc) ~ 10“^ due to finite size effects, the dot¬ 
ted curves do not seem to support power-law scaling of 
This also holds for the broadened ‘box (FFT)’ 
data which establishes that the broadening destroys the 
power law even for uj ^ x- 

To further analyze this surprising finding in a ‘con¬ 
trolled’ setup we took the function f{x) = and 

convoluted it with a Lorentzian of width 10“^. We indeed 
found that for x G [10“^, 10“^] the logarithmic derivative 
of the resulting function shows a behavior quite similar 
to the one of the dotted lines in Fig. [^); in particular, 
it bends away from the expected plateau towards smaller 
values. For larger exponents, let’s say -1/2 as in Figs.|^) 
and b), this does not happen. One is thus tempted to in¬ 
crease the interaction strength and thus 27c —1. However, 
larger interactions also imply larger energy level spacings 
of the (5-peaks and thus require larger x, which cuts off 
the power-law scaling at larger energies. 

As a possible way out we performed a deeonvolution 
of our numerically obtained spectral function along the 
lines of Refs. ED and [32] The resolution was chosen to be 
approximately equal to the broadening %. The deconvo¬ 
lution of numerical data is of course an ill-posed problem 
and the corresponding spectral function can e.g. become 
negative.l^ We were able to perform a stable deconvo¬ 
lution down to {uj — ujmax)/{vFqc) ~ 10“^. In fact, the 
resulting spectra are sujfieiently smooth to perform log- 
arithmie derivatives. These are shown as solid lines in 
Fig. [^). The oscillatory behavior at the lower end sig¬ 
nals the onset of an instability of the deconvolution. The 
deconvoluted data for the box potential lie exactly on top 
of the ones obtained from the unbroadened spectral func¬ 
tion ‘box (ana)’. This gives us confidence that also for the 
other potentials the deconvolution can be trusted. This 
is further supported by the observation that the data for 
the Gaussian potential now clearly support power-law 
scaling with exponent 2xc — 1- For the exponential po¬ 
tential the energy scale cutting off the power-law scaling 
at the ^high energy^ side appears to be of the order 10“^ 
and thus smaller compared to the one of the box (of the 
order 1) and the Gaussian potential (of the order 10“^). 
Therefore no clear plateau is reached for the accessible 
energies. However, the data appear to saturate at the 
expected value 2 jc — 1- 

Along the same lines we next investigate whether or 
not for k — kp > 0 a power law is found close to the 
peak on the inverse photoemission side. In Figs. IT) 
[{k — kF)/qc = 0.01] and c) [{k — kF)/qc = 0.1] we present 
the logarithmic derivative of the deconvoluted spectra to¬ 
gether with ‘box (ana)’ data. The logarithmic derivative 
of the original spectra behave similar to the k = k-p case 
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power law is cut off at {uj — ci;max)/('^F^c) ~ 10 “^ due to 
finite size effects. The deconvoluted box potential data 
fall again exactly on top of the ‘box (ana)’ ones indicating 
that the deconvolution is stable. For (k — /cf)/^c = 0-01 
the data for the Gaussian potential are for sufficiently 
small UJ — cjmax on top of the box potential ones [see 
Fig. [^)]. This does no longer hold for (/c — /cf)/^c = 0.1 
[see Fig. &)] for which the data of the Gaussian poten¬ 
tial are incompatible with power-law scaling. For the 
exponential potential with smaller p this deviation from 
possible power-law behavior sets in at already smaller 
k — k-p [see Fig. |^)]. 

As for the spinful ^ 4 -model the most consistent inter¬ 
pretation of our results is that for any potential with 
p < 00 and fixed k — kp ^ 0 the power-law scaling found 
for the box potential is destroyed by the curvature of the 
potential at g = 0 ; power laws are thus nongeneric. 

We note that this does not only hold when taking cjmax 
as the point of reference. We studied the behavior rela¬ 
tive to other distinguished energies (e.g. Vc[k — k^] and 
the - due to the broadening - apparent threshold). For 
none of these points we find behavior which is consistent 
with power laws. Similarly also the behavior close to 
the threshold on the photoemission side are incompati¬ 
ble with power-law scaling if p < 00 . 

These results imply that the spectral functions result¬ 
ing from the different ad hoc regularizations (not shown 
here), which are all characterized by the threshold power 
laws Eqs. (39) and (40), are nongeneric. 


IV. DISCUSSION 

Our above results for the exact spectral function can 
be summarized as follows: 

1. For k — k-p 0 and generic two-particle potentials 
which are not ‘infinitely fiat’ at momentum transfer 
q = 0 the spectral function does not show power- 
law scaling close to any of the distinguished ener¬ 
gies. Power-law behavior is generically only found 
if all energy scales are sent to zero, e.g. in (/cf, cj) 
for cj ^ 0 (and if ^ 2 ( 0 ) is finite). 


FIG. 4. (Color online) Logarithmic derivative Eq. ( [5^ of 
the spectral function of the spinless g 2 -g 4 -niodel close to the 
inverse photoemission peak for a) (k — kF)/qc = 0 and b) 
{k — kF)/qc = 0.01, and c) {k — kF)/qc = 0.1. The full lines 
labeled by ‘box (ana)’ show the unbroadened results of the 
box potential. The dotted lines in a) result from broadened 
spectra. The other full lines are obtained from the numerical 
spectra after a deconvolution (see the text). The dashed lines 
indicate the expected exponent. The parameters are as in 
Fig.i 


and are thus not shown. As expected the ‘box (ana)’ 
data are consistent with power-law scaling with expo¬ 
nent 7 c — 1 = —7/8 (for the given parameters). The 


2. The ad hoc regularized spectra which are com¬ 
monly studied show finite k — k^ power laws and 
can thus not be considered as generic. The origin 
of this nongeneric behavior is the linearization of 
the spin and charge dispersion relations. The ad 
hoc spectra are plagued by additional artifacts (see 
Sect. [mE ). 

3. The less fiat the potential at g = 0, that is the 
smaller the introduced index p, the faster the dif¬ 
ferences of the spectral line shape compared to the 
one of the box potential - and in many respects 
compared to the one of the ad hoc procedures - 
becomes apparent when k — k-p is increased. Eor 
small p, e.g. the exponential potential with p = 1 , 
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already for (/c — /cF)/gc = 0-1 major differences are 
apparent; see Figs. Ht) and 1^). For /c — /cf 7 ^ 0 
the spectral function is still characterized by the 
dispersing spin and charge peaks. 


To individually study how the curvature of the potential 
at zero momentum transfer modifies the two interaction 
effects of spin-charge separation and power laws with in¬ 
teraction dependent exponents, we have studied the spin¬ 
ful ^ 4 -model and the spinless ^ 2 -^ 4 -model. Combining 
the two limiting cases it is obvious that the same conclu¬ 
sions can be drawn for the full TLM, that is the spinful 
5 ' 2 - 5 ' 4 -model. Explicit results for th e spec tral function of 
the full TLM are presented in Sect. IV C| (see Fig. [^. 

We next discuss the implications of these findings. 


A. Luttinger liquid universality 

1. General eonsiderations 

The TLM forms the low-energy fixed point model un¬ 
der RG flow of a large class of Id correlated fermion 
models.l^ Based on this insight it was suggested that the 
power laws of the ad hoc regularized TLM spectral func¬ 
tion p^{k 7 ^ k-p^) should be observable in other models 
from this class .1^ One can doubt this on general grounds 
as /c — /cf 7 ^ 0 sets a finite scale thus breaking quantum 
critical scale invariance and cutting off the RG flow.l^ Our 
results show explicitly that power laws in {k 7 ^ /cf,cj) 
are indeed not part of LL universality; if they cannot be 
found in the generic TLM, they cannot expected to be a 
universal feature of other models. 

We emphasize that LL universality does not imply that 
for a given model from the LL universality class (other 
than the TLM itself) one simply has to choose proper 
coupling functions ~ which in any case would be 

unknown a priori - and can reproduce details of the spec¬ 
tral function of this model at low energies using the TLM. 
All that is known about the spectral properties of a model 
from the LL universality class is that for u; ^ 0+ the scal¬ 
ing relations p^{k-F^oo) ^ and p^{oj) ^ hold. 

When computing space-time correlation functions of 
the TLM other than the single-particle Green function 
momentum integrals of similar type as encountered here 
appear. To evaluate these ad hoc procedures including 
the linearization of the collective spin and charge disper¬ 
sion are commonly employedOne prominent example 
is the density-density correlation function for momentum 
transfer close to 2 /cfJ^ As for the single-particle spectral 
function after regularization the Fourier integrals can be 
performed analytically leading to power-law scaling of 
the susceptibility with interaction dependent exponents 
close to the characteristic energies :hVc[k^2k-F]. Our con¬ 
siderations imply that also this feature can most likely 
not be considered as being characteristic for LLs in gen¬ 
eral. The density response at vanishing energy and for 
k ± 2 /cf (all energy scales are sent to 0 ) is, however. 


in general characterized by a power-law divergence (for 
repulsive interactions) which indicates the breakdown of 
linear response theory. It is a signature of the sensitivity 
of a LL against single-particle perturbations with mo¬ 
mentum transfer 2 /cf (backscattering) .1^ We note that 
the character of the density-density correlation function 
with small momentum transfer is not modified by the 
momentum dependence of the potentials; it is given by a 
(5-function at energy ujc{q) instead of 


2. Speetra of lattice models 

Directly computing the single-particle spectral func¬ 
tion for translational invariant microscopic lattice models 
of Id correlated fermions still poses a formidable chal¬ 
lenge of quantum many-body theory. Two promising 
routes exists. 

The first one is numerical in nature. Using exact di- 
agonalization (ED), the (dynamical) density-matrix RG 
(DMRG), or different types of quantum Monte-Carlo 
(QMC) approaches, valuable information on the spectral 
function of different models was collected. However, the 
search for power laws requires an exceptional energy res¬ 
olution. This implies that fairly large system sizes and 
low temperatures must be accessible. The model most 
heavily investigated is the Id Hubbard model (and vari¬ 
ants of it). Away from half-filling it is known to fall into 
the LL universality class. ED is restricted to systems 
of a few ten lattice sites which leads to a poor energy 
resolution. Within so-called cluster perturbation theory, 
which is ED-based, it was at least possible to observe 
spin-charge separat ion.l^ The latter was systematically 
studied using QMC.In these studies the finite tem¬ 
perature (and the required analytic continuation) turned 
out to be the main obstacle preventing an analysis of 
p^{k,uj) in terms of power-law scaling. Reference |36] con¬ 
tains the first serious attempt in this direction. The au¬ 
thors use (dynamical) DMRG to obtain broadened spec¬ 
tra. A scaling analysis as a function of the broadening 
was interpreted to be consistent with power-law behav¬ 
ior of the charge and spin peaks. The quality of the data 
is, however, not good enough to either confirm or refute 
power laws at /c — /cf 7 ^ 0 . Generally speaking the nu¬ 
merical results for lattice models from the LL universality 
class are fully consistent with our conclusions. 

A word of warning is in order. As our analysis of 
Sect. |mF] shows it can be very difficult to undoubtedly 
confirm or refute power-law behavior of broadened finite 
size data. This holds even though for the TLM we can 
achieve a comparably high energy resolution and obtain 
data which are sufficiently accurate to employ a very sen¬ 
sitive logarithmic derivative. 

A promising analytical approach to the spectral func¬ 
tion of Id models is built upon the special structure of 
several of the standard lattice models from the LL uni¬ 
versality class, namely the existence of an extensive num¬ 
ber of local integrals of motion. Eor this reason sev- 
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eral models are exactly solvable by Bethe ansatz, which, 
however, does not imply that the single-particle spectral 
function can (easily) be computed exactly.^ In a series 
of papers, see Ref. [37l and references therein, p^{k,uj) of 
the Id Hubbard model was investigated using a ‘pseudo¬ 
fermion dynamical theory’ which is based on the Bethe 
ansatz solution. It was reported that the spectral func¬ 
tion contains power laws even for /c — /cp 7 ^ 0 , with ex¬ 
ponents which depend on k — This finding might be 
related to the ‘nonlinear Luttinger liquid phenomenol¬ 
ogy’ briefly touched below. The spectral function of an¬ 
other exactly solvable Id model, namely the Calogero- 
Sutherland model showing power laws with momentum 
dependent exponents was interpreted in the light of this 
phenomenology.^ 

We emphasize that our findings do not exclude that 
specific models from the LL universality class might show 
finite k — k-p power laws, however, if so for more specific 
reasons than LL universality. Models with equilibrium 
dynamics which are restricted by an extensive number 
of local conserved quantities might be examples for such 
behavior. This is also supported by the observation that 
the power laws reported on in Refs.|T7l|36l and [37l cannot 
only be found at low energies (for very small \k — 
while LL theory is supposed to be applicable in this limit 
only. 


B. Implications for the ‘nonlinear Luttinger liquid 
phenomenology’ 

In an extensive series of papers, reviewed in Ref. [23l 
a novel approach to study corrections to LL behavior 
by the curvature of the single-particle dispersion was de¬ 
veloped. It cumulates in the so-called ‘nonlinear Lut¬ 
tinger liquid phenomenology’. The approach is mainly 
built upon an effective field theory which is motivated 
by lowest order perturbation theory. For the single¬ 
particle spectral function this phenomenology predicts 
(Fermi edge singularity like) power laws with momentum 
dependent exponents even at k — k^ 7 ^ 0. However, the 
field theory relies on the assumption of a momentum in¬ 
dependent bulk interaction (the momentum dependence 
of the interaction with the mobile impurity is consid¬ 
ered) and requires ad hoc ultraviolet regularization. As 
we have shown in (linear) Luttinger liquid theory the 
same assumption leads to power-law behavior which is 
nongeneric rather than universal. This raises the ques¬ 
tion whether the power laws found in ‘nonlinear LL phe¬ 
nomenology’ are robust against a curvature of the bulk 
two-particle potential. 


C. Implications for the interpretation of 
experimental spectra 

When interpreting experimental angular resolved pho- 
toemission data on quasi-Id metallic material^^ the ad 


3 

& 

v + 





w/(uf9c) 



FIG. 5. (Color online) Photoemission part of the spectral 
function of the spinful TLM at {k — /cf)/^c = —O.I. A spin- 
independent interaction with equal inter- and intra-branch 
amplitude is assumed. The interaction strength is chosen 
such that 7c (^ = 0) assumes the values given in the leg¬ 
end. The other parameters are m = 10^ (system size) and 
x/(i’F^c) = 3 • 10“^ (broadening). In a) the spectra for the 
Gaussian potential are shown, in b) the ones for the exponen¬ 
tial potential. The small ‘high-energy’ oscillations for larger 
interactions in a) are a finite size effect. 


hoc regularized spectral function of the TLM is often 
taken paradigmatically. It is e.g. expected that the dis¬ 
tribution of the spectral weight between the spin and 
charge peaks as well as the spacing between the two 
must be exactly as in the analytical expressions given in 
Refs, m [131 HH and |28l A very recent example in which 
this leads to a putative conflict can be found in Ref. |38l 
Our results for the spinful g^- and the spinless ^ 2 -^ 4 - 
models show that the ad hoc regularized spectral func¬ 
tion cannot be considered as universal. Within the TLM 
the details of the distribution of spectral weight, the peak 
distance, and the line shape clearly depend on the form 
of the two-particle potential even when considering small 
Ik — k^l- Such ‘details’ are not part of LL universality and 
expecting quantitative agreement overstresses LL theory. 
We reiterate that all that is known about the T = 0 
spectral properties from the latter is that for cj ^ 0^ the 
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scaling relations p^{kY^oo) ^ and p^{uj) ^ hold. 
The exponent = 7 c + 7 s can be expressed in terms of 
the LL parameters K^/s [see Eq. ([T4|; Kg = 1 for spin- 
rotational invariant systems], whicn also enter in other 
‘observablesIf they can be measured for the same 
system, consistency checks are possible. LL universality 
also makes predictions for the scaling of spectral weight 
as a function T. It was shown that {uj = 0) ^ and 
P>(/cf,cc; = 0) ^ jMEU We note that detail s of the 
analysis of p^{k,uj) at T > 0 for the TLIV#^^ beyond 
the above scaling relation should be taken with caution 
as they rely on an ad hoc regularization procedure. 

It is tempting to employ the (spinful) TLM spectra 
computed with proper foi* comparison to exper¬ 

imental ones beyond the above scaling relations. This 
might lead to a qualitative agreement of certain features. 
One crucial generic feature we found regardless of the q- 
dependence of the two-particle potential are dispersing 
spin and charge ‘peaks’, however, generically not given 
by power-law singularities. To illustrate this we show 
the spectral function of the spinful TLM for the Gaus¬ 
sian and the exponential potential in Figs. Hi) and b), 
respectively. This time we focus on the photoemission 
part p^{k,uj) for {k — kY)/qc = —O.I as this is most easily 
accessible experimentally. The results were obtained by a 
straightforward generalization of the recursive procedure 
put forward for the spinful ^ 4 -model (se e Sect. |IIID I ) 
and the spinless g 2 -g 4 -nciode\ (see Sect. |IIID2p . The 
system size is set by = 10 ^ and the broadening by 


x/('^F^c) = 3 • 10“^. We consider the physically reason¬ 
able case of a spin-independent interaction with an equal 
amplitude of the intra- and inter-branch parts. This im¬ 
plies uJs{(l) = and thus Vs{q = 0) = v-p- The in¬ 

teraction is varied such that 7 c (g = 0) = 1/8, 1/3, and 
9/16. Accordingly, Vc{q = 0 ) increases. As for the above 
discussed limiting cases of the spinful g^- and the spin¬ 
less g 2 -^ 4 -iaodel, the details of the spectral line shape 
obviously depend on the details of the two-particle po¬ 
tential [compare Figs. [^) and b)]. The integral over 
p+(/c,u;) = p^{k,uj) ^{k,uj) at fixed k — k-p must be 
normalized to I which implies that for the exponential 
potential less spectral weight is transferred to p^{k,uj) 
(not shown) than for the Gaussian one. For the former 
the ‘charge peak’ is fairly broad and has an unusal line 
shape. However, the two sets of spectra also share sim¬ 
ilarities. For increasing interaction the spin peak loses 
weight and deforms into a shoulder like structure. To 
summarize, we expect spin-charge separation to be a ro¬ 
bust feature of quasi-Id metalls which should be observ¬ 
able in photoemission experiments. 
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